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We investigate the dynamics of high aspect ratio nanowires trapped axially in a single gradi¬ 
ent force optical tweezers. A power spectrum analysis of the Brownian dynamics reveals a broad 
spectral resonance of the order of a kHz with peak properties that are strongly dependent on the 
input trapping power. Modelling of the dynamical equations of motion of the trapped nanowire 
that incorporate non-conservative effects through asymmetric coupling between translational and 
rotational degrees of freedom provides excellent agreement with the experimental observations. An 
associated observation of persistent cyclical motion around the equilibrium trapping position using 
winding analysis provides further evidence for the influence of non-conservative forces. 


Optical trapping studies of low-symmetry objects, such 
as rods, ellipsoids and oblate spheroids, provides new 
dimensions to the study of optical trapping dynamics. 
In addition to more complex hydrodynamic interactions, 
non-spherical objects are subject to preferential trapping 
orientations and associated restoring torques when dis¬ 
placed away from equilibrium. When the trapping point 
is positioned arbitrarily with respect to the geometric 
centre, perturbative forces can lead to complex cross¬ 
coupling between translational and rotational motion[l|. 
Shape mediated effects have been discussed in the context 
of tailored optical force profiles, including the generation 
of optical lift Q , negative radiation pressure Q , and have 
consequently been implemented in a practical setting for 
passive force clamping in photonic force microscopy [^. 

One interesting peculiarity that has seen growing at¬ 
tention is the tendency of low-symmetry particles to un¬ 
dergo persistent oscillatory motion in thermal equilib¬ 
rium. This type of motion was first articulated for op¬ 
tically trapped rods by Simpson and Hanna, and is re¬ 
lated to the influence of non-conservative optical scat¬ 
tering forces, which induces coupling between rotation 
and translation . For rods aligned nominally along the 
trapping axis (z-axis), the predicted motion is to circu¬ 
late around a central trapping point within a plane (e.g 
XZ plane), coupled with a rocking of the rod orienta¬ 
tion around the orthogonal axis (Oy). The concept has 
its genesis in an earlier study by Roichman et al. who 
explored the influence of the non-conservative scattering 
force in producing circulating currents in the stochastic 
trajectories of spherical particles in an optical tweezers 
0. These concepts have been extended by Saberi and 
Gittes who consider general non-conservative forcing ef¬ 
fects in optical tweezers [3|- 

Our interests lie in experimental investigations of 
weakly tapered, high refractive index, indium phos¬ 
phide (InP) nanowires of uniform composition, held in 
a linearly polarised, gradient force optical tweezers Q. 
Nanowires are a particular subclass of non spherical ob- 
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FIG. 1. (Color online) Schematic of our experimental setup 
and geometry within the optical trap with relevant coordi¬ 
nates and angles. In the experimental arrangement a diode 
pumped solid state laser is used with adaptive optics to create 
a diffraction limited spot within the focal plane of the objec¬ 
tive (OBJ). The nanowires are trapped near the end of the 
wire and aligned along the laser axis. Light scattering from 
the nanowire is collected by the condenser (COND) and used 
to perform back focal plane interferometry with a position 
sensitive detector (PSD) 


jects with nanoscale radial dimensions and lengths that 
extend to several micrometers, giving aspect ratios in 
excess of 100:1. For trapping experiments, the extended 
axial dimension can exceed the depth of focus of the trap¬ 
ping beams by several times the Rayleigh range, thereby 
creating a trapping force that is localised to a particular 
point on the axis. Earlier studies have shown that the 
equilibrium trapping position of these nanowires has been 
determined to be near one end, with the major segment 
of the length sitting below the trap far from the geo¬ 
metric centre [^, which is qualitatively consistent with 
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modelling of tapered rods 0 - The optically trapped 
nanowires are aligned closely to the optical axis of the 
tweezers and experience small perturbations from the 
equilibrium point due to stochastic motion. Note that 
these trapping conditions exclude the possibility of ma¬ 
jor orientation changes that are observed in other studies 

Eiiia 

Here we present experimental evidence for the presence 
of resonant fluctuations in the stochastic motion of opti¬ 
cally trapped nanowires. A power spectral density anal¬ 
ysis of the nanowire motion reveals a broad resonance 
peak of the of the order of a kHz, which is strongly de¬ 
pendent on trapping power. Modelling of the dynamical 
equations of motion of the trapped nanowire with asym¬ 
metric cross coupling terms between translational and ro¬ 
tational degrees of freedom provide good agreement with 
the experimental observations, suggesting that the effect 
is due to the action of non-conservative forces. Further, 
a winding number analysis of the Brownian reveals clear 
cyclical motion around the trapping point, which is an¬ 
other characteristic property of non-conservative effects. 

Our experimental apparatus consists of a single beam 
gradient force optical trap generated at the focus of a 
high numerical aperture oil immersion microscope ob¬ 
jective lens (Nikon E Plan Achromat lOOx 1.25NA) us¬ 
ing a Nd:YAG laser emitting at 1064nm (Laser Quan¬ 
tum IR Ventus). Dynamic steering of the laser focus 
is provided by a two-axis acousto-optic deflector (AOD) 
(Gooch and Housego 45035 AOBD) and spatial light 
modulator (SLM) (Hamamatsu LGOS-SLM xl0468-03) 
positioned within the optical train at planes conjugate 
with the back aperture of the objective. SLM-based aber¬ 
ration correction and overfilling of the back aperture are 
used to produce a diffraction limited spot optimised for 
optical trapping. The power of the trapping laser at the 
imaging plane is controlled by the AOD and set between 
10 mW and 90 mW. A schematic diagram of the trapping 
arrangement is shown in FIG [T] 

The indium phosphide (InP) nanowires used in the 
study are epitaxially grown on InP (lll)B substrates by 
metal-organic chemical vapour deposition (MOCVD) us¬ 
ing gold nanoparticle catalysts. The nanowires are be¬ 
tween 5 and 10p,m have a tapered aspect of approxi¬ 
mately 0.3° as characterized by electron microscopy. Fur¬ 
ther information on the growth process can be found else¬ 
where, [13 1 . 

For tweezers measurements, the nanowires are dis¬ 
persed into solution by sonication and transferred to a 
hermetically sealed chamber at a dilution where single 
nanowires can be isolated. Visualising the nanowires 
in the trap is achieved using a charge-coupled device 
(GGD) camera (AVT Stingray) together with diascopic 
illumination supplied by the condenser. Back focal 
plane interferometry is performed by collecting the for¬ 
ward scattered light from the trapped nanowires with 
a high numerical aperture condenser lens (Olympus S 


Plan Fluor 40x 0.6NA) and projected onto a position 
sensitive diode (PSD) (Pacific Sensor, DL16-7PGBA) lo¬ 
cated at a conjugate plane to the condenser back aper¬ 
ture. The PSD is connected to a field-programmable gate 
array (FPGA) multi-function acquisition card (National 
Instruments PGIe-7852R) and time series measurements 
of the nanowire position fluctuation with respect to the 
trapping centre are recorded at a sampling rate of lOkHz 
for a period of 10 seconds. The linearity of the PSD de¬ 
tector response with respect to the nanowire position in 
this orientation has been previously verified [ 8 |. 
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FIG. 2. (Color online) (a) Power spectral density of a single 
trapped nanowire exhibiting resonant oscillations (red) and 
the fitted spectrum (black) accounting for rotation / trans¬ 
lational coupling and associated non-conservative effects, (b) 
Power spectrum of the same nanowire trapped close to the 
cover slip (blue) where the resonance significantly suppressed. 
The dotted (red) line is a Lorentzian fit and the solid (black) 
is the model, (c) Autocorrelation function of the nanowire 
motion at different heights showing resonant (red) and non¬ 
resonant behaviour (blue). The inset is a bright field image 
of the un-trapped nanowire. 


The stochastic motion of a nanowire within an op¬ 
tical trap is normally described by a set of uncoupled 
Langevin equations in terms of translational (x, y, z) and 
angular coordinates {9x,0y). Each coordinate has an as¬ 
sociated linear restoring torque/force, drag force and in¬ 
dependent noise terms, Qit) with zero mean and variance 
{C,i{t)C,i(t -b t)) = 2kbTjiri(T). The trapping experiments 
are performed in a low Reynolds Number aqueous envi¬ 
ronments and the associated drag coefficient, 7 , is given 
by, 
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Here 7 _l and 7 || are the drag coefficients perpendicular 
and parallel to the symmetry axis of the cylinder, and 
je is the rotational drag coefficient, ijo is the viscosity of 
water dj_, (5|| and Sg are correction factors accounting for 
the ends of the cylinder which depends on the ratio of the 
length (/) to diameter (2r) of the cylinder [lj-[l6|. Power 
spectrum analysis of this system yields a characteristic 
Lorentzian dependence for each degree of freedom where 
the corner frequency is given by the fc = Kij2'K^i , where 
Ki is the optical spring constant. 

A typical power spectrum of a trapped nanowire (blue) 
is presented in FIG[2Kb), together with the corresponding 
Lorentzian fit (dotted). Using the physical dimensions of 
50 nm by 5 /rm in length, (FIG [2Kc) [inset]), we estimated 
Kx to be 2.9 pN/im“^mW“^ in this measurement, which 
is comparable to previously published results. Solving 
the Langevin equation for rotation yields a similar form of 
the power spectrum from which the restoring torque con¬ 
stant can be extracted. From previous work analysing au¬ 
tocorrelation data it is noted that the characteristic fre¬ 
quency for rotation is somewhat smaller compared with 
the translation motion and thus is not resolvable in these 
measurements 
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When the trapping height above the cover slip is in¬ 
creased a broad resonant peak emerges in the power spec¬ 
trum (red), shown in FIGUJa). The resonance peak per¬ 
sists only as the nanowire is trapped at a critical distance 
of about 50 /im from the cover slip, and reverts back to 
the normal response with decreased trapping height. The 
emergence of the resonance is also reflected in the au¬ 
tocorrelation function of the nanowire trajectory, where 
a fluctuation between positive and negative correlation 
can be seen, as shown in FIG [SJc). The full height of 
the chamber is estimated to be in excess of lOO/xm and 
so we exclude the possibility of direct interaction with 
the chamber walls. We note that changing the height 
of the trap will induce spherical aberrations and accen¬ 
tuate other aberrations already present within the trap. 
Increased aberrations will change the relative contribu¬ 
tions from gradient and scattering forces for a given fixed 
input power, and we hypothesise that this is the origin 
of the height dependence of the resonance. 

In order to account for the emergence of this resonant 
peak we consider a more general form of the equations 
of motion which include cross-coupling terms. The dy¬ 
namical equations in an over-damped environment can be 
written as a single 5-coordinate vector equation in terms 


of the generalized forces (neglecting rotation around the 
z-axis), 


-T^+at)-Kq = Q. ( 2 ) 

Here, q = {x,0y,y,Ox, z) are the generalised co¬ 
ordinates, F is the diagonal hydrodynamic drag matrix, 
C (t) is the Langevin force and K is the trap stiffness ma¬ 
trix. 
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Importantly, we only include off-diagonal terms in the 
stiffness matrix that couple between translation and ro¬ 
tation within a plane, i.e. a rotation in the xy plane {6y) 
will induce a translation in x. This is justified as the non¬ 
conservative forces (radiation pressure) that are driving 
the coupling are principally directed along the z-axis and 
thus should induce rotation around the orthogonal axis. 
The above-mentioned approximations leads to two cou¬ 
pled equations of motion that retain the non-conservative 
effects, and that can be treated independently. An anal¬ 
ysis of the two-coordinate stochastic equations leads to a 
power spectral density for each coupled coordinate (e.g. 
i = x,9y) that has the form 18|: 


or ^ VrUJ + aj , . 

27r((w2 + y) 

where fj, and e are the respective trace and determi¬ 
nant of the matrix M = r~^K, and rji is associated with 
the noise strength. The term ax = g + 

is related to the cross-coupling matrix elements in M. 
In interpreting the solution in terms of the experimen¬ 
tal system we note that if the coupling relates to the 
non-conservative scattering force, the off-axis compo¬ 
nents of the stiffness matrix should be non-symmetric, 
e.g. Kg^^x ^ ^x,ey The form of the power spectral 
density is similar to that of a damped harmonic oscilla¬ 
tor, such as is observed in aerosol tweezing experiments 
[l9| . This is to be expected as both represent coupled 
first order differential equations; in the case of the un¬ 
derdamped oscillator coupling is between restoring and 
inertial terms, whilst in the present case it is between 
translation and rotation. A fit of the power spectrum 
presenting the resonant peak in FIG [IJa) using Eq. |4] 
shows excellent agreement with the experimental data. 
In fitting the data we apply an additional offset to ac¬ 
count for detector noise. Indeed a fit of the data taken at 
the lower trapping height FIG [2Kb) also reveals a strongly 
suppressed peak at a similar frequency. 
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On inspecting the fitting parameters between of the 
different trapping heights we make the important empir¬ 
ical observation that a non-zero resonance peak is only 
observed under conditions where the one of the cross¬ 
coupling terms is negative. Such a scenario is ensured in 
this system as the trapping point is far from the centre of 
rotation - the action of rotation under Brownian motion 
will induces a translation at the trapping point which is 
in the opposite direction to the associated restoring force. 



FIG. 3. (Color online) Power spectral density of a single 
nanowires held at a fixed trapping power measured for differ¬ 
ent trapping powers. The peak frequency shifts given by the 
increases linearly with trapping power, whilst the Brown¬ 
ian fluctuations are strongly suppressed. The solid lines rep¬ 
resent the fitted model which includes non-conservative cross¬ 
coupling between translation and rotation. 

To further investigate the dynamics of the resonance, 
the power spectrum of a trapped nanowire of length 
4.6/xm was recorded as a function of input trapping laser 
powers ranging from 10 mW to 90 mW. For low symme¬ 
try particles the trap stiffness scales linearly with power. 
The power dependent results are presented in FIG [3] 
along with the associated fits using Eq. H] again we note 
exceptional agreements between experimental results and 
the model. Qualitatively we see that the peak position 
of the resonance is blue shifted from 200 Hz to 1800 Hz 
with increasing trapping power together with a reduction 
of the amplitude of the fluctuations. The parameters as¬ 
sociated with the fits are given in FIGlU We observe that 
e increases quadratically with trapping power, whilst /i 
increases linearly. As the trap stiffness is proportional 
to trapping power and /r = tr{M) is proportional to the 
trap stiffness we should expect a linear increase in this 
parameter. Similarly, the as e is the determinant of M 
and should therefore be proportional to the square of the 
matrix elements, and hence proportional to the square 
of the power, ij is roughly constant at across all trap¬ 
ping powers, which indicates that there is no obvious 
heating across the different trapping powers. Finally, we 


note that whilst the trapping power dictates the relative 
strengths of the trap stiffness (conservative) and radi¬ 
ation pressure (non-conservative), it should not modify 
the ratio of the two competing forces under conditions 
where the harmonic approximation for the trap stiffness 
is valid. 




FIG. 4. (Golor online) Fitting parameters of the power spec¬ 
tral density data presented in FIG O (a) e increase quadrat¬ 
ically and /i linearly with trapping power, which corresponds 
to linear increase in both the peak frequency and width, (b) 
The noise strength, rj, shows no clear dependence on trapping 
power indicating a constant temperature, whilst a increases 
quadratic with trapping power. All dependencies are consis¬ 
tent with a trap stiffness that is linearly proportional to the 
trapping intensity in the presence of asymmetric coupling. 

In additional to the characteristic power spectrum de¬ 
scribed above, another important property of this system 
is its penchants to produce a continuous cycling of the po¬ 
sition/orientation under steady state conditions. In order 
to discern the presence of cyclic motion a winding num¬ 
ber analysis of the nanowire trajectories is performed. 
In our case the winding events are denoted when the 
nanowire trajectory moves around the centre of the trap 
in either a clockwise (negative) or anticlockwise (positive) 
motion. In an unbiased system the rates of cycling in ei¬ 
ther direction is balance leading to an average of zero 
winding. The results presented for different powers in 
FIG[5Ka), indicate that nanowires do indeed have a ten¬ 
dency to persistently cycle around a central point within 
the trap in a particular sense; short timescale fluctuations 
correspond to Brownian fluctuations in the trend. The 
winding behaviour is persistent over many seconds and 
is reproducible for each individually trapped wire. The 
sense of the cycling tends to be in a particular sense, sug¬ 
gesting an underlying bias within the trapping geometry, 
however winding in the opposite sense is also observed 
for some trapping events. We find that for different pow¬ 
ers the rate of winding increases before plateauing at 
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higher powers; this is clearly observed in a plot of the 
gradient of the winding number with trapping power, 
shown in FIG El^b). The winding property is dissipa¬ 
tive in nature and is sustained only by the action of the 
non-conservative scattering forces. We may understand 
the process most intuitively by considering the following: 
when the nanowire is tilted with respect to the direction 
of the scattering force (along the z-axis) it will experi¬ 
ence a greater pushing force. If the nanowire position 
lies below the focus it will be pushed through the focus 
to a position that is balanced by the axial restoring force. 
The additional lift force will only be alleviated when the 
tilt is reduced leading to a rotation back to the equilib¬ 
rium position. The nanowire is then free to drop below 
the equilibrium position along the z-axis. 
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FIG. 5. (Color online) (a) A winding number plot of 
the nanowire motion for different trapping powers. Small 
timescale fluctuations in the winding decorate a clear mona- 
tonic increase with time, indicating persistant cycling in one 
direction, (b) The rate of winding, given by the gradient 
of the winding plots, is observed to increase approximately 
trapping power for lower powers and plateau at higher pow¬ 
ers. The winding also takes on a distinctly nonlinear trend at 
very high trapping frequencies. 

In conclusion, we provide clear experimental and theo¬ 
retical evidence for the influence of non-conservative cou¬ 
pling between translation and rotational modes in opti¬ 
cally trapped nanowires. We see experimentally that this 
leads to a distinct resonance peak in the power spectrum 
and an accompanying winding of the particle trajecto¬ 


ries. This model system provides exciting new oppor¬ 
tunities to study complex dynamic behaviour involving 
non-conservative forces and may lead to exciting new in¬ 
sights into analogous biophysical systems. 
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